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Q-i, In this letter the explicit form of evolution operator of three atoms Tavis- 



Cummings model is given, which is a generalization of the paper quant-ph/0403008. 

US 

The purpose of this letter is to give an explicit form to the evolution operator of Tavis- 
Cummings model (PP) with some atoms. This model is a very important one in Quantum 
Optics and has been studied widely, see [2] as general textbooks in quantum optics. 

We are studying a quantum computation and therefore want to study the model from 
this point of view, namely the quantum computation based on atoms of laser-cooled and 
trapped linearly in a cavity. We must in this model construct a controlled NOT gate or 



* E-mail address 
t E-mail address 
* E-mail address 
§ E-mail address 
^E-mail address 



fujii@yokohama-cu.ac.jp 

s035577d@yokohama-cu.ac.jp 

s035559g@yokohama-cu.ac.jp 

suzukita@gm.math.waseda.ac.jp 

s035588a@yokohama-cu.ac.jp 



1 



other controlled unitary gates to perform a quantum computation, see [3] as a general 
introduction to this subject. 

For that aim we need the explicit form of evolution operator of the models with one, 
two and three atoms (at least). As to the model of one atom or two atoms it is more or 
less known (see j3]), while as to the case of three atoms it has not been given as far as we 
know. Since we succeeded in finding the explicit form for three atoms case we report it . 

The Tavis-Cummings model (with n-atoms) that we will treat in this paper can be 
written as follows (we set h = 1 for simplicity). 

An n 

H = uj\ l a* a + - a ? ] ® 1 + W +) ® a + a t ] ® o f ) , (1) 
^ i=i i=i 

where uo is the frequency of radiation field, A the energy difference of two level atoms, 
a and a) are annihilation and creation operators of the field, and g a coupling constant, 
and L = 2 n . Here a\ + \ a\ ' and are given as 

(j\ s) = 1 2 ® - • • ® 1 2 ® a s <S> h ® ■ • • ® h (* - position) e M(L, C) (2) 

where s is +, — and 3 respectively and 
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Here let us rewrite the hamiltonian (0) . If we set 

Th Tt "1 Tl 

s+ = E^ (+) , s- = £*{-\ s 3 = ~$>l 3) , (4) 

j=l i=l i=l 

then can be written as 

H = ool L ® a f a + A£ 3 ® 1 + g (S+ ® a + S_ ® a 1 ") = F + V, (5) 



which is very clear. We note that {S + , S~, S 3 } satisfy the su(2)-relation 
[S 3 ,S + ]=S + , [S 3 , SJ\ = -S-, [S+, = 2S 3 . (6) 

1 J.C.Retamal et al might have obtained the same result by another method |S] 



However, the representation p defined by p(u+) = S + , p(cr_) = p(a 3 /2) = S 3 is a 
reducible representation of su{2). 

We would like to solve the Schrodinger equation 

ijU = HU=(H + V)U, (7) 

where U is a unitary operator (called the evolution operator). We can solve this equation 
by using the method of constant variation. The result is well-known to be 



U(t) = (e~ itujSa (g) Q~ itLua ' <a ^j Q -itg{s+®a+S-®a<) 



under the resonance condition A = ou, where we have dropped the constant unitary 
operator for simplicity. Therefore we have only to calculate the term (JBJ) explicitly, which 
is however a very hard task 2 . In the following we set 

A = S + <g) a + S- <8) a ] (9) 

for simplicity. We can determine e~ lt9A for n = 1 (one atom case), n = 2 (two atoms 
case) and n = 3 (three atoms case) completely. 

One Atom Case In this case A in Q is written as 

/ 



A 1 

By making use of the relation 



a 
at 



(10) 



A, 2 



' aa) \ f N+l 
v a) a ) \ N 



11) 



with the number operator iV we have 

/c OS ( i9 VivTT) -i!=&gE 



(12) 



2 the situation is very similar to that of the paper quant-ph/03 12060 in [Jj 
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We obtained the explicit form of solution. However, this form is more or less well-known, 
see for example the second book in [2] . We note that (fT2j) can be decomposed as 



cos (tgy/N) 
( cos (tgy/N + 1 



1 



. tan(tgVN) t 







V 








1 



COS^tg^) J 



.tan^tgy/W+i) 



(13) 



This is a Gauss decomposition of unitary operator. This may be used to construct 
a theory of "quantum" representation of a non-commutative group, which is now under 
consideration. 



Two Atoms Case In this case A in 



|9j) is written as 

a a » 

at a 

a f a 

at a) / 



(14) 



Our method is to reduce the 4 x 4-matrix A2 in (JT4J) to a 3 x 3-matrix Si in the 
following to make our calculation easier. For that aim we prepare the following matrix 
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then it is easy to see 





y/2a! 
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where B\ = J + <g> a + J_ ® a* and { J + , J_} are just generators of (spin one) irreducible 
representation of Q- We note that this means a well-known decomposition of spin 



I®! 



0©1. 



Therefore to calculate e lt9A2 we have only to do e ttgBl . Noting the relation 



/ 



Bi 



2(2iV + 3) 



2(2JV + 1) 



we obtain 



-itgB\ 



2{2N - 1) i 
\ 



B = DBi, 



hi bi 3 

&21 ^22 ^23 
^ &31 ^32 ^33 ) 



(15) 



where 



6n 

&13 



iV + 2 + (JV + l)cos (tgy/2(2N + 3) 

2iV + 3 
-1 + cos (tgy/2(2N + 3) 
2ivT3 



&12 = -i- 



sin W2(2iV + 3) 



V2F+3 



a 2 , &21 = — 



sin (tgy/2(2N + l)) f 
V2A^+ 1 



6 22 = cos W2(2JV+1) , 6 23 = -? 



sin (W2(2JV+1) 



'31 



'33 



-1 + cos (tgyj2(2N-V 
2N-1 



(a 1 ") 2 , 6 32 = -i 



V2iV+T 

sin (tgy/2(2N- 1 



y/2N-\ 



N - 1 + iVcos (tgyj2(2N - 1) 



2iV- 1 



5 



Three Atoms Case In this case A in @ is written as 



/ 



A, 



a a 
at 
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at at a 
atOOOOaaO 
at at a 
at a f a 
at at at 



(16) 



We would like to look for the explicit form of solution like (fT2*|) or ()15j) . If we set 
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then it is not difficult to see 

^ a 
at 



T ] AvT 
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at 





V^at 
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\/3at 
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This means a decomposition of spin | Cg> | £g> | = \ © \ © §. Therefore we have only to 
calculate e~** sB2 , which is however not easy. The result is 

-itgB2 



e 

^ /i(JV + 2) -V3i/i 2 (^V + 2)a 2^/3^ + 2)a 2 -6i/i 3 (iV + 2)a 3 ^ 



■y/3ih 2 (N + l)a* / 2 (W+1) -2^ 3 (^ + l)a 2^/3(^ + 1; 



a 



2 



2v / 3/ 3 (iV)(a t ) 2 -2i/i 3 (^V)a t / 4 (JV) -y/Zih A {N)a 

y -Qth 3 (N - l)(at)3 2v/3/ 3 (iV-l)(at) 2 -y/3ih 4 (N - l) a t f 5 (N-l) ) 



where 



/i(JV) = |v+(iV)co8(^^/A + (iV)) - V _(iV)cos(^^/A_(iV))} /(2y/d(N)), 
f 2 (N) = \w + (N)cos(tgy/X + (N)) - w.(N)cos(tgy/X_(N))\ /(2y/d(N)) 



h 2 (N) = \^^sm(tg^N)) - ^^ S m(tg^X~{N)) ) /(2^d(N)), 
[y/X + (N) J\-(N) J 

f 3 (N) = (cos(tg^X + (N)) - coBitgy/X-iN))} /(2y/d(N)), 



= - 7 =L=sin(^ V / A + (7V)) - - 7 =L=sin(^ V / A_(iV)) /(2y/d(N)), 

/ 4 (JV) = {mA0cos(^v/A-(A0) " u-Wcos^^/MiV))} /(2\AW) ? 

MAO = (^S=sin(^ V / Al(A 7 )) - -^^sin(^^V))} /(2y/d{N)), 
[y/X-(N) v/A + (iV) J 

/ 5 (7V) = {w + (N)cos(tgJxZ(N)) - W _(W)cos(t<^A + (W))} /(2y/d(N)), 
h 3 (N) = y\ + (N)sm(tgy/\ + (N)) - Jx^)sm(tgJx~(N))} /(2y/d(N)) 



and 



(17) 



d(N) = 16N 2 + 9, A± = 5iV ± Jd(N),v± = -2N - 3 ± \MA0, iu± = 2iV - 3 ± Jd(iV). 



We obtained the explicit form of evolution operator of the Tavis-Cummings model 
for three atoms case, so there are many applications to quantum optics or mathematical 
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physics, see for example 4J. In the near future we will apply the result to a quantum 
computation based on atoms of laser-cooled and trapped linearly in a cavity [H] . 

We conclude this paper by making a comment. The Tavis-Cummings model is based 
on (only) two energy levels of atoms. However, an atom has in general infinitely many 
energy levels, so it is natural to use this possibility. We are also studying a quantum 
computation based on multi-level systems of atoms (a qudit theory) Therefore we 
would like to extend the Tavis-Cummings model based on two-levels to a model based 
on multi-levels. This is a very challenging task. 
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